Numerical black holes: A moving grid approach.

Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein’s equations with the harmonic slicing condition is presented. Significant errors in the mass function are shown to arise from the steep gradient zone behind the black hole horizon, which challenge the Computational Fluid Dynamics numerical methods used in the code. The formalism is extended to moving numerical grids, which are adapted to follow horizon motion. The black hole exterior region can then be modeled with higher accuracy. PACS numbers: 04.20.Cv Typeset using REVTEX On leave of absence from Departament de F́ısica. Universitat de les Illes Balears. E-07071 Palma de Mallorca, SPAIN.